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Abstract. The angular momentum projected Generator Coordinate Method, with the quadrupole moment 
as collective coordinate and the Gogny force (DIS) as the effective interaction, is used to describe the 
properties of the ground state and low-lying excited states of the even-even Neon isotopes ^'^"^'*Ne, that is, 
from the stability valley up to the drip-line. It is found that the ground state of the N=20 nucleus ^"^Ne is 
deformed but to a lesser extent than the N=20 isotope of the Magnesium. In the calculations, the isotope 
^^Ne is at the drip-line in good agreement with other theoretical predictions. On the other hand, rather 
good agreement with experimental data for many observables is obtained. 
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1 Introduction 



■ The properties of the ground and lowest-lying excited states 
' of nuclei close to the stability valley are determined to a 

' great extent by the underlying mean field in which all the 

■ nucleons move. The nuclear mean field, obtained through 
[ the mean field approximation to the nuclear many body 
, problem, provides the important concepts of magic num- 
' bers (or shell closures) and that of spontaneous symmetry 
] breaking. For nuclei with proton and/or neutron numbers 

■ close to the magic ones symmetry conserving (i.e. non su- 
perconducting and spherical) ground states are expected. 
On the other hand, nuclei away from the magic configu- 
rations are expected to show strong symmetry breaking 
at the mean field level leading to deformed (and super- 
conducting) ground states which are the heads of bands 
generated kinematically by restoring the broken symme- 
tries -like the well known example of the rotational bands. 

Almost thirty years ago, experimental evidences were 
found in the neutron-rich light nuclei around '^^Na p]|2j 
pointing to the breaking of the N=20 magic number. The 
first attempt to understand the phenomenon from a theo- 
retical point of view p] suggested that correlations beyond 
mean field coming from the restoration of the rotational 
symmetry could be strong enough as to overcome the 
mean field shell effects and lead to deformed ground states 
for closed shell nuclei Hke '^^Na and '^^Mg. Further theo- 
retical studies using the Shell Model (SM) approach were 
carried out fl|0 and it was suggested that an intruder 
configuration consisting of a two particle two hole neu- 
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tron excitation from the sd shell into the fy /2 one was the 
responsible for deformation in the ground state of ^^Mg. 

The recent availability of Radioactive Ion Beam fa- 
cilities in several laboratories like Ganil, GSI, MSU and 
Riken to cite a few and the development of very efficient 
mass separators and solid state detectors has made possi- 
ble to measure up many properties concerning the ground 
state and the lowest lying excited states of many exotic, 
neutron-rich light nuclei. In particular, the exploration of 
both the = 20 and TV = 28 shell closures far from 
stability has proven to be a rich source of new phenom- 
ena. Among the variety of available experimental data, 
the most convincing evidence for a deformed ground state 
in the region around = 20 is found in the "^^Mg nu- 
cleus where both the excitation energies of the lowest lying 
2+0 and 4+ states and the B{E2, 0+ 2+) transition 
probability have been measured. The low excitation en- 
ergy of the 2+ state, the high value of the B(E2) transi- 
tion probability and also the ratio i?(4+)/i?(2^) = 2.6 
are fairly compatible with the expectations for a rota- 
tional band. Additional evidence comes from the neigh- 
boring isotope ^''Mg where the £^(2^^) is only 0.66 MeV, 
the B{E2,0+ 2+) is 631 e^fm^ and the E{At)/E{2t) 
ratio is 3.2, that is, the 0+, 2^ and 4+ satisfy all the re- 
quirements to belong to a strongly deformed rotational 
band ptTHj . 

From a theoretical point of view, the ground state of 
■^^Mg is spherical at the mean field level. However, when 
the zero point rotational energy correction (ZPRE) stem- 
ming from the restoration of the rotational symmetry is 
considered, the energy landscape as a function of the quad- 
rupole moment changes dramatically and ■^^Mg becomes 
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deformed [3lllllll2lll3l[Ti] . A more careful analysis of the 
energy landscape including the ZPRE correction reveals 
that, in fact, there are two coexistent configurations (pro- 
late and oblate) with similar energies indicating thereby 
that configuration mixing of states with different quad- 
rupole intrinsic deformation has to be considered. There- 
fore, an Angular Momentum Projected Generator Coordi- 
nate Method (AMP-GCM) calculation with the quadru- 
pole moment as collective coordinate is called for. We have 
appHed this method in Ref. jH] to the study of the nuclei 
30-34-|y|-g ^j^jj ^jjg Gogny force. We have obtained prolate 
ground states for ^^"^'^Mg indicating that the TV = 20 
shell closure is not preserved for the Mg isotopes. More- 
over, a good agreement with the experimental data for 
the 2+ excitation energies and B{E2) transition probabil- 
ities was obtained. The method has been used to study 
quadrupole collectivity of the Si isotopes ^21) the N=28 
isotopes as well as the superdeformed band in ^^S 
with quite a good success. 

The purpose of this paper is to extend the AMP-GCM 
calculations to the study of the Neon (two proton less than 
Magnesium) isotopes from A=20 up to A=34. Contrary 
to the Magnesium isotopic chain, the Neon neutron drip 
line (N=22) is very close to the neutron magic number 
N=20 and therefore the study of the Neon isotopes is a 
good testing ground to examine both the systematic of de- 
formation and the possible erosion of the N=20 spherical 
shell closure very close to the neutron drip fine. Another 
interesting point is the study of the possible magicity of 
N=16 suggested in recent analyses [19] . 

The paper is organized as follows: In section 2 a brief 
overview of the theoretical framework is presented. In sec- 
tion 3.1 the mean field results are discussed. In the next 
section the effect of angular momentum projection on the 
mean field observables is described. Finally, in section 3.3 
the results of the configuration mixing calculations are 
presented and compared to the experimental data and 
other theoretical approaches. We end up with the con- 
clusions in section 4. 



2 Theoretical framework 

As mentioned in the introduction ours is a mean field 
based procedure where the underlying mean field is deter- 
mined first and then additional correlations beyond mean 
field are included. Those additional correlations are han- 
dled in the framework of the angular momentum pro- 
jected Generator Coordinate Method (AMP-GCM) with 
the mass quadrupole moment as generating coordinate. As 
we restrict ourselves to axially symmetric configurations, 
we use the following ansatz for the X = AMP-GCM 
wave functions 

K) = J dq20fUl20)P^0 l^feo)) . (1) 

For each angular momentum / the different AMP-GCM 
states (labelled by a) are linear combinations of the set 
of angular momentum projected intrinsic wave functions 



\f{(l2a)) generated by solving the Hartree-Fock-Bogoliubov 
(HFB) equation constrained to yield the desired mass quad- 
rupole moment 520 = {vi<l2o)\ - l/2(a;2 + y^) |(y9(g2o))- 

The intrinsic wave functions are restricted to be axially 
symmetric (i.e. if = 0) and are obtained by solving the 
HFB equation with the Gogny interaction [201 (DIS pa- 
rameterization [2T|). The HFB equation is discretized by 
expanding the quasiparticle operators associated to the 
intrinsic wave functions |<y3((j'2o)) in a Harmonic Oscillator 
(HO) basis containing eleven major shells and with equal 
oscillator lengths (in order to make the basis closed under 
rotations [22]). As we are deahng with quite fight systems 
we have to consider the center of mass problem. This is 
handled by subtracting the center of mass kinetic energy 
both in the calculation of the energy and in the HFB vari- 
ational procedure. Finally, concerning the Coulomb inter- 
action, we have only taken into account its contribution to 
the direct field in the variational procedure. The exchange 
Coulomb energy (computed in the Slater approximation) 
is added, in a perturbative fashion, at the end of the cal- 
culation and the contribution of the Coulomb interaction 
to the pairing energy is completely disregarded. 

In order to obtain the angular momentum projected 
wave functions we use the standard angular momentum 
projector operator restricted to iiT = states [221 

n'o = ^^^^ / dnSMe-^^^' e-/34e-.7i. . (2) 

Finally, the "collective amplitudes" fl{q2o) as well as the 
energies of the AMP-GCM states \<Pi) are obtained through 
the solution of the Hill- Wheeler (HW) equation 

I dq',on'{q20,q'2o)fU'l2o)^Ei J dq^o (920 , 920 ) /J (920 ) 

(3) 

which is given in terms of the projected norm 

A^'('?20,92o) = (^(920)1 P^o 1^(920)) (4) 
and the projected hamiltonian kernel 

(920,920) = (^(92o)|i?Po'o 1^(920)) • (5) 

As the generating states Pqq \'f{q2a)) are not orthogonal, 
the "collective amplitudes" /i(92o) cannot be interpreted 
as probability amplitudes. Instead, one usually introduce 
[2i| the "collective" wave functions 

5^(920) - J dq',of Uq2o) (A^' (920, 9^0)*)'^' (6) 

which are orthonormal / d92offi*(92o).9a'(92o) Sjj^S^^^' 
and therefore their module squared has the meaning of a 
probability. 

In order to readjust the particle number on the average 
for the projected wave functions the Hamiltonian in Eq. 
© has been replaced by H' = H-Xn{N-No)-Xp{Z-Zo) 
where A„ and Ap are chemical potential parameters (see 
Ref [25, and references therein). 
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Fig. 1. In panels a) and d) the ground state HFB energies 

of 32,34 

are plotted as functions of the HO length 6o (&x = hz = bo) 
for the Nsheii ~ 8, 9, 10, 11, 12, 13 and 17 bases. The curve corresponding to the basis used in the present work {Nsheii ~ 11) is 
plotted as a dashed line. In panels b) and e) the HFB energies of ^^•^''Ne computed with bo = 2.1fm and Nshdi = ll(full line) 
and Nsheii = 12 (dashed line) are plotted as functions of the quadrupole moment. In panels c) and f), the energy differences 
between Nsheii = H and Nsheii = 12 calculations are plotted as functions of the quadrupole moment. 



The B{E2) transition probabilities are computed using 
the AMP-GCM wave functions 

B(£;2,/,-.//) = ^^x (7) 
I dq2odq'2ofaf*iq2o){IfQ'20 W Q2 \ \ ^i92o)/i: (^20) 

with 

(If 120 \\Q2 \\I^q20) ^^fh 2 // \ 

(2/, + l)(2// + l) Z^y^fifiO) 

r d(3 sin f3d'l^,{P){^{q'2„) I Q2,e~'^'- \ ^(<72o)) 
Jo 

where the indices i and f stand for the initial and final 
states and 02^^ are the charge quadrupole moment op- 
erators. As we are using the full configuration space no 
effective charges are needed. Further details on the com- 
putational procedure can be found in Ref |25| . 



3 Discussion of the results 
3.1 Mean field approximation. 

Before discussing the mean field results, the convergence of 
our calculations with the size of the Harmonic Oscillator 
(HO) basis used to discretize the HFB equation has to 
be tested as we are dealing with near drip-line nuclei like 
32,34^g_ First, one should note that, due to the proximity 
of the drip-line, the full HFB approximation must be used 
[2filt27| . It is also evident that absolute convergence for the 
binding energies can only be attained for HO basis with an 
infinity number of shells {Ngheii)- However, it is expected 
that other physical observables Hke excitation energies, 
transition probabilities, etc can be accurately described 
with a finite number of shells. 

In order to determine how many shells are needed for 
an accurate description of the physical observables we have 
first studied the behavior of the HFB energy as a function 
of the oscillator length parameter bo {b± = bz — bo) and 
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the number of shells for the nuclei '^^Ne and Ne. The re- 
sults are plotted in panels a) and d) of Fig.Q] As expected 
the curves become more and more flat for increasing val- 
ues of N shell and already the N shell = 17 basis can be 
considered as a good approximation for an inflnite basis 
in both isotopes as the dependence of the energies on the 
oscillator length is very weak for a wide range of 60 val- 
ues. However, Nsheii = 17 makes the HO basis too big 
for our purposes as we have to evaluate the AMP-GCM 
hamiltonian kernel which is a quantity requiring two or- 
ders of magnitude more computing time than the whole 
HFB calculation. We have found that N shell = 11 is a 
good compromise between accuracy and computing time 
as the energy still shows a weak dependency on the os- 
cillator length around the minima located at 5o = 2.1 fm 
for both '^^Ne (panel a) and '^^Ne (panel d). As an exam- 
ple of the adequacy of the N shell = 11 calculations let us 
consider the two neutron separation energy of ■^"'Ne. The 
Nsheii — 17 ground state energies for '^^Ne and ^''Ne are 
403 keV and 420 keV lower than the N shell = 11 energies 
implying that the two neutron separation energy of ■^'^Ne 
computed with N shell = 11 differs by 17 keV from the one 
computed with Nsheii = 17. 

In order to study the suitability of the Nsheii = 11 
basis for other nuclear properties we have performed cal- 
culations with Nsheii = 12. The mean fleld results for both 
calculations are compared on the right hand side panels of 
Fig.^ In panels b) and e) the energy landscapes of '^^Ne 
and '^'^Ne are shown as a function of the quadrupole mo- 
ment for the Nsheii = 11 and Nsheii — 12 calculations. In 
the upper panels c) and f) we have represented the energy 
differences between both calculations. The conclusion is 
that in the region —0.7b < 520 < l-2b the shape of the 
energy landscapes do not change much when the basis is 
increased from Nsheii = H to Nsheii = 12. As we will see 
later on this is the range of 920 values where the collective 
dynamics is concentrated and therefore it is not expected 
to find significant differences between the Nsheii = 11 and 
Nsheii = 12 results. 

In Fig.[2the mean field potential energy surfaces (MF- 
PES) are plotted as a function of the axially symmetric 
quadrupole moment (720 for the even-even Neon isotopes 
20-34]v^g The MFPES shown do not include the Coulomb 
exchange energy and they have been shifted to accommo- 
date them in a single plot (see figure caption). 

Both ^°Ne and ^^Ne are prolate deformed in their ground 
states. In ^''Ne the prolate ground state corresponds to 
920 = 0.4b (/?2 = 0.37) and an oblate local minimum also 
appears at 920 = —0.1b {(32 = —0.09) with an excitation 
energy of 1.71 MeV. In the case of ^^Ne the ground state 
corresponds to 920 = 0.5b {(32 — 0.40) and another local 
minimum is found at (720 = —0.2b {P2 = —0.17) with an 
excitation energy of 2.24 MeV. The nucleus ^^Ne is a clear 
example, in the considered isotopic chain, of very strong 
shape coexistence since, while the oblate ground state is 
locate at q2o = —0.3b {P2 — —0.22), a prolate isomeric 
state is also found at (720 = 0.2b {P2 = 0.15) with an ex- 
citation energy with respect to the oblate ground state 
of 77 keV. On the other hand, the nuclei ^^"^°Ne show 
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Fig. 2. Mean field potential energy surfaces for the considered 
Neon isotopes, plotted as a function of the axially symmetric 
quadrupole moment. The curves have been shifted to show 
them in a single plot. The corresponding energy shifts are given 
in the plot. 
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Fig. 3. In each panel, proton (thick full lines) and neutron 
(thick dashed lines) particle-particle are de- 

picted as a function of the quadrupole moment. Also the energy 
differences £'jjp]3(92o) — -Eg.s., where Eg.s. corresponds to the 
ground state HFB configuration, are plotted as thin full lines. 
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spherical ground states indicating that the N=20 shell clo- 
sure is preserved at the mean field level. The MFPES of 
both ^^'^^Ne are particularly flat around their spherical 
ground states. In the nucleus ^°Ne we obtain a prolate 
shoulder at 920 = 0.8b {P2 = 0.37) at an excitation energy 
of 2.70 MeV with respect to the spherical ground state. 
This prolate shoulder is around 1 MeV higher than the 
one found in similar HFB calculations in'^^Mg (see for ex- 
ample EHj). In the drip-line systems '^^Ne and '^^Ne, pro- 
late deformed ground states are found. The ground states 
have 920 = 0.6b {P2 = 0.26) and gao = 0.8b {p2 = 0.31), 
respectively. In addition, an oblate isomeric state is found 
in •^''Ne at 920 = —0.3b {(32 = —0.12) with an excitation 
energy of 2.39 MeV with respect to the prolate ground 
state. 

Another interesting point concerns the stability of the 
Ne isotopes against neutron emission. From the absolute 
mmimum of the MFPES' depicted in Fig. [H (the one of 
^°Ne is marked with a horizontal line) it can be deduced 
that '^^Ne is not stable against two neutron emission. This 
is in clear contradiction with the experimental results as 
they indicate that '^^Ne is indeed stable against two neu- 
tron emission [28,29,30]. On the other hand, the nucleus 
^''Ne is slightly stable against two neutron emission, but 
it is not against four neutron emission. 

In Fig. the proton and neutron particle-particle cor- 
relation energies —Epp — ^Tr (Ak*), which are commonly 
used to discuss pairing correlations (see for example Refs. 
|2()II31) 1. are plotted as a function of the quadrupole de- 
formation for all the isotopes considered. The evolution of 
the particle-particle correlation energies is well correlated 
with the structures found in the MFPES. Non-zero proton 
pairing correlations are found in all the spherical or oblate 
minima. In addition, sizeable neutron pairing correlations 
are found in ^"'^^Ne and ■^^''^^Ne for the spherical and the 
oblate minima. Vanishing proton pairing correlations are 
found in the prolate side in a window starting at 0.5b and 
ending in 1.0 b in ^°Ne. For the other isotopes the start- 
ing and ending points increase with the neutron number. 
Neutron pairing correlations vanish in the ground states 
of both ^^Ne and '^°Ne. In the later case this is a conse- 
quence of the N=20 shell closure found at the mean field 
level. 

It should be stressed here that the unphysical collapse 
of pairing correlations found in Fig.|3|is an indication that 
one should also consider in these isotopes dynamical pair- 
ing correlations and their coupling to the quadrupole de- 
gree of freedom in the scope of a formalism beyond the 
mean field in order to treat them in an equal footing. 
However, this multidimensional configuration mixing cal- 
culation is a cumbersome task with the present-day com- 
putational facilities. 



3.2 Correlations beyond the mean field: Angular 
momentum projection. 

Before considering the full AMP-GCM it is instructive to 
look into the angular momentum projected energy sur- 




-1.0 0.0 1.0 2.0 

£720 (b) 



-1.0 0.0 1.0 2.0 

£720 (b) 



Fig. 4. Angular momentum projected potential energy sur- 
faces (full lines) for the nuclei ^'^~"^'*Ne and for the spin values 
I'^ — 0"'", 4"^and 6"^, plotted as a function of the axially 
symmetric quadrupole moment q2o- The mean field potential 
energy surfaces are also plotted as lines with boxes. In each 
nucleus, the energies are referred to the energy of the — 0"*" 
ground state. For the nuclei ^^"^"^Ne we have also included 
(dashed lines) the projected results corresponding to the cal- 
culation with Nsheii = 12. 



faces (AMPPES) defined as 

^'(920) 



'^^(920,920) 
A/"^ (920, 920) 



(8) 



and shown in Fig. 31 for the nuclei ^'^^^''Ne and / = 0, 2, 
4, and 6. The corresponding mean field energy landscapes 
(dashed curves with boxes) are also included for compari- 
son. For details on the missing points in the 7 = 2,4 and 6 
curves refer to (2^. The most remarkable fact about Fig. 
H is how strongly the restoration of the rotational sym- 
metry modifies the mean field picture of the 7 = con- 
figurations. Contrary to the mean field case, two minima, 
one prolate and the other oblate, are obtained for all the 
Neon isotopes for = 0+. The prolate minimum is, with 
the exception of ^''Ne, the absolute minimum in all the 
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Fig. 5. A comparison of the angular momentum projected 
energies of the '^"Ne and '^^Mg nuclei as functions of the quad- 
rupole moment for = C^and 2"^ (referred to the — 
0"*" absolute minimum in each nucleus). 



isotopes considered. For increasing spin values either the 
energy difference between the prolate and oblate minima 
increases or the oblate minimum is vi^ashed out. 

The nucleus ^^Ne, with its I'^ — 0+ oblate ground 
state is the exception in all the Neon isotopes studied. 
The orbital responsible for such an oblate minimum is 
the neutron d| orbital which becomes fully occupied in 
this nucleus and favors oblate deformations. The absolute 
minimum remains oblate for = 2+ but already at I"^ — 
4+ it becomes prolate deformed. 




0.0 1.0 2.0 -1.0 0.0 1.0 

£720 (b) £720 (b) 

Fig. 6. Band diagrams for the nuclei ^""^''Ne obtained by 
solving the Hill- Wheeler equation of the AMPGCM. The 7 = 
AMPPES are also plotted to guide the eye. See text for further 
details. 



3.3 Correlations beyond the mean field: Angular 
momentum projection and configuration mixing. 



In addition, shape coexistence is expected in the nuclei 
26Ne, 2*Ne and 3°Ne as their r = 0+ prolate and oblate 
minima are very close in energy (414, 577 and 936 keV, re- 
spectively). These minima are separated by barriers which 
are 2.4, 2.2 and 1.3 MeV high, respectively. 

Finally, it is also worth to compare the intrinsic quad- 
rupole deformation of the ground state of '^°Ne with the 
one of ^^Mg In Fig. El we have plotted the r = 0+ 
and 2+ AMPPES for both ^ONe and ^^Mg. From this plot 
we conclude that the absolute minimum of the — O"*" 
AMPPES in 3°Ne has half the deformation of the cor- 
responding minimum in ^^Mg. On the other hand, the 

= 2+ deformations are practically identical. From this 
results one may expect, as is the case, very similar spec- 
troscopic quadrupole moments for the first 2+ states in 
^"Ne and '^^Mg and a reduction of the B(E2) transition 
probability in '^"Ne as compared to the one of ^^Mg. 



The AMPPES of the previous subsection show for some 
nuclei and/or some spin values the phenomenon of shape 
coexistence and therefore configuration mixing has to be 
considered in order to gain a better understanding of the 
structure of these states. We have considered configuration 
mixing in the framework of the Angular Momentum Pro- 
jected Generator Coordinate Method (AMP-GCM) de- 
scribed in section |2l The intrinsic axial quadrupole mo- 
ment (720 with values in the range —1.5b < (720 < 2.5b 
and with a mesh size Aq2o of lOfm^ has been chosen as 
generating coordinate. 

In Fig. El we have plotted the AMP-GCM energies E'^ 
obtained by solving the Hill- Wheeler equation and placed 
them along the (720-axis according to their "average intrin- 
sic quadrupole moment" defined as 

920^ / ^920 I 9iiq2o) P 920 (9) 
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Fig. 7. Collective wave functions squared for the ground states (a = 1) and the spin values r = 0+, 2+, 4+and 6+ of the nuclei 
20-34]\jg 'pijg corresponding projected energy curve is also plotted for each spin value. The y— axis scales are in energy units 
and always span an energy interval of 13 MeV (minor ticks are 0.5 MeV apart). The collective wave functions have also been 
plotted against the energy scale after proper scaling and shifting, that is, the quantity + 15x | gi{q2o) \^ is the one actually 
plotted. With this choice of the scales we can read from the figure the energy gain due to quadrupole fluctuations by considering 
the position of the wave functions' tail relative to the projected curve. 



In this figure we have also plotted the r = 0+ AMPPES 
to guide the eye. The first noticeable fact is that the cor- 
related 0+ ground states have a lower energy than the ab- 
solute minimum of the corresponding AMPPES. This en- 
ergy gain due to the quadrupole correlations increases the 
binding energies and can be relevant for a proper descrip- 
tion of the two neutron separation energies. The second 
prominent feature is that configuration mixing decreases 
the deformation of the 0^ ground states with respect to 
the minimum of the AMPPES. The ground states of the 
nuclei ^'*Ne and ^^Ne become spherical whereas the ones 
of ^^Ne and ■^"Ne become weakly deformed. The other 
nuclei remain well deformed in their ground states and 
develop a rather well defined rotational band up to the 
maximum spin considered for tr = 1. In addition, a well 
defined rotational band is obtained for / > 2 in "^^Ne. 
On the other hand, the excited states (cr — 2) only show 
a rotational band pattern for those nuclei well deformed 
in their ground state. Another interesting point concerns 
the energy gain due to configuration mixing of the ground 
state: it is 0.6 MeV for ^°Ne, increases up to 0.9 MeV for 
^''Ne and ^^Ne and then monotonically decreases up to 
0.55 MeV in '^'^Ne. The energy gain is correlated with the 
shape of the AMPPES for 1=0, namely, if the AMPPES 
shows a well defined minimum the energy gain is smaller 
than in the case where the minimum is broader. 

To gain a deeper insight into the intrinsic structure 
of the AMP-GCM sates the ground band (cr = 1) col- 



lective wave functions squared | gi{q2o) P are plotted in 
Fig. □ along with the corresponding AMPPES. For ^°Ne 
and ^^Ne the ground state collective wave functions are 
well inside the prolate wells indicating that both systems 
are dominated by prolate deformations in the considered 
spin range. The average deformation values are 0.39b, 
0.58b, 0.59b and 0.58b for the 0^", 2+, 4+ and 6^" states 
in ^°Ne, while the corresponding values in ^^Ne are 0.39b, 
0.56b, 0.61b and 0.63b. In the nucleus ^''Ne, the configu- 
ration mixing calculation provides a spherical 0^ instead 
of the oblate absolute minimum found in its AMPPES. 
An almost spherical 2f state (920'^''^"^ = —0.08b) is also 
found in this nucleus. On the other hand, a band cross- 
ing takes place for the state 4^ and the collective wave 
function becomes prolate deformed with gjcT — 0.39b 
and ^2(7^ '^"^ = 0.79b. Similar results have recently been 
found ,25.J in the iV = 14 system ^^Mg. The experimental 
spectroscopic quadrupole moment of the lowest 2+ state 
in ^®Mg is -13.5 efm^ indicating that this is a pro- 
late state. However, the low excitation energy of the Oj 
in ^^Mg (3.588 MeV) compared with the same quantity in 
^■^Mg (6.432 MeV) clearly indicates strong shape coexis- 
tence between oblate and prolate solutions. Unlike ^^Mg, 
there is not experimental information concerning the spec- 
troscopic quadrupole moment of the lowest 2+ state in 
^^Ne but again the experimental excitation energies of 
the Oj in 24Ne and ^'^Ne (4.764 MeV and 3.691 MeV, 
respectively) are significantly lower than the 6.235 MeV 
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Table 1. Ground state spectroscopic quadrupole moments 
Qspec^j^^ = 1) in eW for r = 2+, 4+ and 6+ in the nu- 
clei ^'^~^'*Ne. Experimental data |82IM4) are shown in boldface 
whereas Shell Model results |35I36| are shown in brackets. For 
details see the main text. 
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Fig. 8. The two neutron separation energies obtained with the 
different approximations considered in this paper are compared 
to experimental values taken from Ref. (37ii38i| . and the QMCD 
results of Ref. |39I40| 



measured in ^^Ne [32ll33| pointing towards a strong shape 
coexistence in the ground state of the former nuclei. Our 
results predict a strong shape coexistence for the 0^ and 
2f states in ^"'Ne as well as for the states in ^^Ne and 
^®Ne that manifest itself in the excitation energies of 
3.509 and 3.005 MeV for ^4.26^5^ respectively. 

The 0^ wave functions in both ^^Ne and ^*Ne show a 
great admixture of prolate and oblate configurations that 
leads to spherical ground states on the average. On the 
other hand, for > 2+ the collective wave functions in 
these two nuclei become prolate deformed with deforma- 
tion values of 0.34b, 0.59b and 0.94b for the states 2-^, 4f 
and 6^ in ^^Ne while the corresponding values in ^®Ne are 
0.26b, 0.93b and 1.04b. Our results do not fully support 
the suggestion of ^Hj concerning the magicity of N=16 
in neutron rich light nuclei. It is true that ^^Ne has a 
spherical ground state, but the deformation of the 2^ is 
too strong as to be considered a vibrational state. Similar 
results are also found in Ref. (2^] for the N=16 nucleus 

28Mg. 

The 0^ wave function in ^°Ne also shows a significant 
admixture of the oblate and prolate configurations and, 
as a consequence, the deformation in the ground state is 
reduced to 0.16 b that represents one third of the value 
corresponding to the absolute minimum in the I'^ = 0"*" 
AMPPES. This clearly shows that, as in ^^Mg pltT5ll25l . 
the deformation effects in ^°Ne are the result of a sub- 
tle balance between the zero point corrections associated 
with the restoration of the rotational symmetry and the 
fiuctuations in the collective parameters (in our case the 
axially symmetric quadrupole moment). Prom the com- 
parison with the value of g2(7" "^"^ = 0-43 b already found 
[T5II25] in ^^Mg we conclude that dynamical deformation 
effects are strongly suppressed in ^°Ne as could have been 
forecasted from the different AMPPES topology we have 
already seen in FigEl On the other hand, the 2^, 4^ and 
6f wave functions in '^°Ne are inside the prolate wells and 
the average deformations of 0.83b, 0.93b and 0.94b are 
very close to the ones found in '^^Mg. 



In both Ne and Ne, the ground state collective 
wave functions become prolate. The dynamical deforma- 
tion values for the O;,^, 2^, 4^" and 6f in ^^Ne are 0.42b, 
0.90b, 0.98b, 1.05b while the corresponding values in '^^Ne 
are 0.69b, 0.94b, 0.99b and 1.0b. All these values show 
the stability of deformation effects in Neon isotopes as we 
move towards the drip-line. 

Here we will make a few comments on the AMP-GCM 
results in the nuclei '^^Ne and '^'*Ne when the basis is in- 
creased from N Shell = 11 to N shell = 12. The effect on 
the projected energy landscapes can be seen in the cor- 
responding panels of Fig. 21 where Nsheii = 12 curves are 
plotted as dashed lines. The results with Nsheii = 12 and 
N Shell = 11 only show very small differences at large abso- 
lute values of 520- On the other hand, in the region of phys- 
ical significance (—.7b < (720 < 1.2b) the results are prac- 
tically indistinguishable. Since this range of (720 values is 
the one where the collective dynamics is concentrated (i.e, 
the tails of the collective wave functions go to zero outside 
this range) we do not expect big changes in the excitation 
energies. For the ■^^Ne nucleus, it turns out that the ex- 
citation energies of all the members of the ground state 
rotational band obtained in the Nsheii = 12 calculation 
are around 7 keV higher than the ones in the Nsheii = 11 
case, i.e., only the 0^ state has been pushed down. For 
•^^Ne the average shift is 6 keV. As a consequence, the 
transition gamma ray energies remain unaltered by the 
increase of the basis size. On the other hand, the excita- 
tion energies (with respect to the true ground state) of 
the members of the excited rotational band decrease on 
the average 40 keV in ^^Ne and 32 keV in ^''Ne and there- 
fore, as in the previous case, the intra-band gamma ray 
energies remain the same. From these results we conclude 
that our calculations are well converged in terms of the 
basis size. 

Coming back to the discussion of our results, in Ta- 
ble m we present our results for the ground band spec- 
troscopic quadrupole moments. A very good agreement is 
observed between the calculated spectroscopic values of 
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Fig. 9. On the left panel the experimental values |38l82l41l42l48l44j of the B(£2,0+ 2+) transition probabilities are 
compared to our predictions and other theoretical results |39l4fl| . On the right panel, the same comparison is made for the 2^*" 
and 4^ excitation energies. 



the 2J" states in ^°Ne and ^^Ne and the experimental val- 
ues (shown in boldface) for these nuclei which are -23 ehx? 

and -17 efm^ j^^l, respectively. In addition, our results 
are consistent with the Shell Model predictions (shown in 
brackets) for the spectroscopic quadrupole moments of the 
2^ states [^Kli.Sfij . Obviously, the discrepancies are larger 
for the shape coexistent nuclei ^'*~^*Ne, It is also worth 
mentioning here that the value obtained for the spectro- 
scopic quadrupole moment of the 2^ state in ■^"Ne (-15.80 
efm^) is only slightly smaller than the similar quantity in 
^^Mg (-19.50 efm^). 

In Fig.lHlwe compare the results of the AMP-GCM two 
neutron separation energies S2N = {N — 2) — {N) 

with the corresponding mean field results (see subsection 
13.111 and also with the available experimental values taken 
from Ref. [S3ISHI- The AMP-GCM improves the S2N of 
^^Ne as compared to the mean field results and also makes 
bound (unHke the mean field prediction) the nucleus ■^^Ne 
in good agreement with the experimental results f2RI29l 
I30j ■ On the other hand, the nucleus '^^Ne becomes unstable 
against two neutron emission in the AMP-GCM approach. 
Our results are very similar to the ones predicted in the 
framework of the QMCD 39,40j specially around TV = 20 
but we differ in the prediction concerning ■^'^Ne. Finally, it 
is worth mentioning that the AMP-GCM binding energy 
is the sum of the mean field binding energy of the intrin- 
sic state plus the energy gain due to the restoration of 
the rotational symmetry plus the energy gain due to the 
configuration mixing. Therefore, the differences in S2N ob- 
tained in the AMP-GCM and the mean field are due to 
the last two contributions. The analysis of those contri- 
butions shows that the rotational energy correction is the 
main responsible for the differences observed in S2N and 
therefore is the ingredient needed to make '^^Ne stable. 



In Fig.inithe excitation energies of the 2^ and 4^ states 
and the B{E2,0f — > 2^) transition probabilities obtained 
in the AMP-GCM are compared with the available exper- 
imental values and also with the predictions of the QMCD 
!39"40]. Concerning the B{E2,0t 2f) transition prob- 
abilities we clearly see, from panel a), that the agreement 
with the available experimental data |33ll32ll41ll42ll43ll44j 
is rather satisfactory and in most cases (with the excep- 
tion of ^^Ne where our prediction appears a little bit low) 
our results stay within the experimental error bars. On 
the other hand, our results are also consistent with the 
predictions of the QMCD [211101 • Our results although 
not as good as the QMCD ones, are very satisfactory con- 
sidering that the parameters of the Gogny force have not 
been fitted to the region and/or the physics of quadrupole 
collectivity and also that no effective charges have been 
used in our calculations of the transition probabilities. 

Our results for the 2f excitation energies agree well 
with the experiment in ^"Ne and ^^Ne |33ll32] and also 
with the theoretical result of Otsuka in ^''Ne. For the 
other isotopes our 2+ excitation energies are always higher 
than experiment and the QMCD predictions. The same 
happens for the 4+ excitation energies. The disagreement 
between our results and the experiment (or the QMCD 
predictions) can be attributed to the fact that the nu- 
clei involved are nice examples of shape coexistence. In 
the presence of shape coexistence the quadrupole degree 
of freedom probably is not enough to describe accurately 
the observables and other degrees of freedom like triaxi- 
ality or pairing fiuctuations should be included. However, 
and based on the nice agreement between our electromag- 
netic transition probabilities and the experimental ones, 
we can conclude that the quadrupole degree of freedom is 
the main ingredient in the physical picture of the neutron 
rich Neon isotopes. 



10 



R.R. Rodriguez-Guzman, J.L. Egido and L.M. Robledo: Quadrupole collectivity of neutron-rich Neon isotopes. 



4 Conclusions 

We have performed angular momentum projected Genera- 
tor Coordinate Method calculations with the Gogny inter- 
action DIS and the mass quadrupole moment as generat- 
ing coordinate in order to describe quadrupole collectivity 
in the even-even nuclei ^°~^^Ne. The lighter isotopes ^°Ne 
and ^^Ne as well as the heavier ones ^^Ne and ^^Ne are well 
deformed in their ground states whereas the other isotopes 
are either spherical (^''Ne and ^^Ne) or slightly deformed 
like ^*Ne and ^"^Ne. The later isotope is less deformed in 
its ground state than its isotone -^^Mg. The two neutron 
separation energies compare well with experimental data 
and show that '^^Ne is bound. Moreover, the B{E2) tran- 
sition probabilities from the ground state to the 2+ state 
are well reproduced in all the cases. Only the 2+ excitation 
energies come too high as compared to the experiment for 
the shape coexistent isotopes. 

This work has been supported in part by the DGI, Ministerio 
de Ciencia y Tecnologia (Spain) under project BFM2001-0184. 
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